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Abstract. Multiplicity distributions of hadrons produced in central nucleus-nucleus 
collisions are studied within the hadron-rcsonancc gas model in the large volume limit. 
In the canonical ensemble conservation of three charges (baryon number, electric 
charge, and strangeness) is enforced. In addition, in the micro-canonical ensemble 
energy conservation is included. An analytical method is used to account for resonance 
decays. Multiplicity distributions and scaled variances for negatively charged hadrons 
are presented along the chemical freeze-out line of central Pb+Pb (Au+Au) collisions 
from SIS to LHC energies. Predictions obtained within different statistical ensembles 
are compared with preliminary NA49 experimental results on central Pb+Pb collisions 
in the SPS energy range. The measured fluctuations are significantly narrower than a 
Poisson reference distribution, and clearly favor expectations for the micro-canonical 
ensemble. 
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1. Introduction 



For more than 50 years statistical models of strong interactions [H [2, [3] have served 
as an important tool to investigate high energy nuclear collisions. The main subject 
of the past study has been the mean multiplicity of produced hadrons (see e.g. Refs. 
[U El El [£])• Only recently first measurements of fluctuations of particle multiplicity 
[8] and transverse momenta [9] were performed. The growing interest in the study 
of fluctuations in strong interactions is motivated by expectations of anomalies in 
the vicinity of the onset of deconfinement [TO] and in the case when the expanding 
system goes through the transition line between quark-gluon plasma and hadron gas 
[TT] . In particular, a critical point of strongly interacting matter may be signaled by 
a characteristic power- law pattern in fluctuations [12]. Apart from being an important 
tool in an effort to study the critical behavior, the study of fluctuations within the 
statistical hadronization model constitutes an essential test of its validity. 

Fluctuations are quantified by the ratio of variance of a multiplicity distribution 
to its mean value, the so-called scaled variance. There is a qualitative difference in 
the properties of mean multiplicity and scaled variance of multiplicity distributions in 
statistical models. In the case of mean multiplicity results obtained within the grand 
canonical ensemble (GCE), canonical ensemble (CE), and micro-canonical ensemble 
(MCE) approach each other in the large volume limit. One refers here to as the 
thermodynamical equivalence of statistical ensembles. It was recently pointed out [13J, 
that corresponding results for the scaled variance are different in different ensembles, 
and thus the scaled variance is sensitive to conservation laws obeyed by a statistical 
system. The differences are preserved in the thermodynamic limit. 

In this contribution we briefly summarize recent results [H] on multiplicty 
fluctuations in the hadron resonance gas model. 

We will first discuss a simple example, accessible to analytical solutions, and then 
sketch how to generalize the procedure to a general multi-specie hadron gas. Lastly a 
comparison of model calculations to recently released NA49 data [15J on charged particle 
multiplicity fluctuations is shown. 



2. Ultra-relativistic Gas of Neutral Particles 



The GCE partition function of an ultra-relativistic gas composed of only neutral 
Boltzmann particles is given by: 
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where V is the volume of the system, T its temperature, and g the degeneracy factor 
due to the particles internal spin. The number of GCE micro-states with fixed particle 
number N is given by the Fourier integral over the generalized GCE partition function: 
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where the Wick-rotated fugacity e l ^ N is introduced to fix particle number N. The 
normalized multiplicity distribution, i.e. the probability to find the system in a state 
with exactly N particles is then given by the ratio: 

z n (v,t) (vgi) N ( T 3 ^ 
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Similarly one can find the number of micro-states in GCE with fixed particle number 
N and fixed energy E: 
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where Z MCE (V, E, N) denotes the MCE partition function. The number of 
micro-states with fixed energy E, but arbitrary particle number iV is simply 
Z E {V,T) = £ Z E ' N (V, T) in GCE, or Z MCE (V, E) = E Z M (V, E, N) in MCE. 

N=l N=l 

The difference is just that the former includes the Boltzmann weight e~ B//T , i.e. 
Z E (V,T) = Z MCE (V,E) e~%. Consequently the MCE multiplicity distribution is 
given by [161 E] : 
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Phrased differently, the Pmce(N) can be expressed as the conditional distribution 
Pgce(N\E) = P GCE (N, E)/P GCE {E). 

The scaled variance of a multiplicity distribution is defined as: 
(N 2 ) - (N) 2 
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It can be shown that in the large volume limit both distributions ([3]) and (jSJ) have the 
same mean value (N) , but rather different scaled variance [16] . The scaled variance in 
MCE converges to ujmce — 0.25, while in GCE it is equal to one, ojgce = 1 (Poisson 
distribution). 



3. Generalization to Multi-Specie Ideal Gas 

Above procedure can be easily generalized to a multi-specie hadron resonance gas. When 
enforcing conservation of three Abelian charges, Qi = (B,S,Q), as well as energy and 
three momentum, E k = (E, P x , P y , P z ), the MCE partition function is given by 
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where the single particle partition function is: 
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Particle number 'conservation' would result in an additional integral. 

Above relation has two implications. The first one is technical. In MCE calculations 
one has to deal with a heavily oscillating (or even irregular) integrand. The integrants 
of Eqs. (|H ED are however very smooth. For large volume the main contribution comes 
from small region around the origin, and an analytical expansion is therefore possible 
US]. Z® J > E \V,T, fij) converges then to a Multivariate-Normal-Distribution, while finite 
volume corrections can be given in the form of Hermite polynomials of low order |17j . 
The analytical approximation is valid only around the maximum of the distribution. 

The second implication is more subtle. Rather than calculating the proper MCE 
partition functions, one can also, as in Eq.Q, define Pmce(N) through a joint 
distribution of energy, momentum, charges and particle number in GCE. Pmce{N) 
would then be called the conditional probability distribution, while Pgce(N) is called 
the marginal distribution of the partition function (JTj). 

4. Hadron-Resonance Gas Model and NA49 data 

In order to predict the energy depedence of P(N) along the chemical freeze-out line 
for Pb+Pb (Au+Au) collisions, the freeze-out parameters should be expressed as a 
function of collision energy. Here we follow the procedure described in [UJ [20]. The 



dependence of hb on the collision energy is parameterized as |6j: /ig (v^vtv 



1.308 GeV • (1 + 0.273 y / ijv r /v) _1 , where the cm. nucleon-nucleon collision energy, 
y / sjv r /v, is taken in units of GeV . The system is assumed to be net strangeness 
free, 5 = 0, and to have the charge to baryon ratio of the initial colliding nuclei, 
Q/B = 0.4. These two conditions define the strange, {is, and electric, /iq, chemical 
potentials. Temperature is defined by the condition: average energy per particle is 
equal to 1 GeV [21]. Finally, the strangeness saturation factor, 75, is parameterized 
as [7] 75 = 1 — 0.396 exp (— 1.23 T/fis). This determines all parameters of the 
model. Although various parameterizations are known in the literature [51 E], resulting 
differences for multiplicity fluctuations are small [14J. 

For the calculations, an extended version of the THERMUS framework [22] was 
used. The THERMUS particle table includes all known hadrons and resonances up 
to a mass of about 2.5 GeV and their respective decay channels. We use quantum 
statistics, but disregard Breit-Wigner width of resonances. We assume that in the 
studied reactions the system volume is large and finite volume corrections to fluctuations 
and mean values can be neglected. Resonance decay has been dealt with as described 
in [T71 [20] . The effect of finite experimental acceptance was taken into account by 
an 'uncorrelated particle' approximation [131 EH 120]. The results for multiplicity 
fluctuations of negatively charged hadrons in the full momentum space and a comparison 
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Figure 1. Left: The scaled variances for negatively charged particles, ui~ , both 
primordial and final, along the chemical freeze-out line for central Pb+Pb (Au+Au) 
collisions. Different lines present the GCE, CE, and MCE results. Symbols at the 
lines for final particles correspond to specific collision energies ranging from SIS to 
LHC energies. The arrows show the effect of resonance decay. Right: The lines 
show acceptance corrected values for lu~ in the SPS energy range. The points show 
preliminary data of NA49 [T5]. Total (statistical +systematic) errors are indicated. 
Figures are taken from [14] 
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Figure 2. The ratio of the multiplicity distributions to a Poisson reference distribution 
with same mean value for negatively charged hadrons produced in central (1%) Pb+Pb 
collisions at 20 A GeV, 30 A GeV, 40 A GeV, 80 A GeV, and 158,4 GeV (from left to 
right) in the NA49 acceptance [T5]. Preliminary experimental data (solid points) of 
NA49 [T5] are compared with prediction of the hadron-resonance gas model obtained 
within different statistical ensembles, the GCE (dotted lines), the CE (dashed-dotted 
lines), and the MCE (solid lines). Figures are taken from [T3] 



to data in the NA49 acceptance are shown in Fig. [IJ The multiplicity distributions in 
the NA49 acceptance for central Pb + Pb collisions at 20, 30, 40, 80, and 158AGeV 
together with the predictions for the GCE, CE, and MCE versions of the model are 
shown in Fig. [2l 

The measured multiplicity distributions are significantly narrower than the Poisson 
one and allow to distinguish between model results derived within different statistical 
ensembles. The data agree surprisingly well with the expectations for the micro- 
canonical ensemble and exclude the canonical and grand-canonical ensembles. 
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